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Polymer  electrolyte  membrane  fuel  cells  are  proving  to  be  a  clean  and  efficient  source  of  energy.  Nowa¬ 
days,  extensive  research  efforts  are  being  focused  on  bringing  this  technology  to  everyday  use.  An 
important  aspect  when  integrating  fuel  cells  in  practical  applications  is  their  ability  to  respond  to  load 
demand.  With  respect  to  this,  due  to  their  complex  internal  dynamics,  fuel  cells  belong  to  the  group 
of  more  slowly  responding  sources.  In  order  to  make  them  more  generally  applicative  they  are  often 
connected  with  a  battery  or  a  super-capacitor  via  a  power  converter  to  form  a  hybrid  power  source.  A 
control  algorithm,  designed  for  such  a  system,  represents  an  interesting  challenge:  it  has  to  adapt  to 
varying  working  conditions  and  operate  optimally  in  terms  of  efficiency  and  reliability,  while  minimiz¬ 
ing  any  impacts  on  the  degradation  of  the  components.  Here,  we  present  an  approach  using  supervisory 
control  automaton  that  switches  between  the  system’s  operational  modes  and  sets  the  references  for  the 
lower-level  control  loops.  The  evaluation  of  the  efficiency  and  degradation  is  carried  out  in  a  simulation 
using  a  model  of  the  widely  used  1.2-kW  Ballard  Nexa  power  module. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  with  a  membrane  made  from  a  polymer  material  (Poly¬ 
mer  Electrolyte  Membrane  fuel  cells,  PEM  FC)  are  a  clean  and  efficient 
source  of  electrical  energy  [1-3].  However,  despite  the  first  practi¬ 
cal  implementations  dating  back  to  the  era  of  manned  flights  into 
space  in  the  early  1960s,  they  have  only  recently  become  the  sub¬ 
ject  of  intensive  study  and  application;  in  particular  in  the  fields  of 
mobile  applications,  uninterruptible  power-supply  systems,  auxil¬ 
iary  power  units  and  vehicle  drives.  The  main  advantages  of  fuel 
cells  are  that  the  energy  storage  is  in  the  form  of  hydrogen  (so 
the  energy  is  not  lost  over  time),  a  high  efficiency,  a  low  operating 
temperature  and  a  high  power  density  [4,5]. 

In  a  FC  stack  there  is  an  ongoing  exogenous  electrochem¬ 
ical  reaction,  where  hydrogen  and  oxygen  (from  the  air)  are 
transformed  into  water;  at  the  same  time  energy  is  released  in  the 
form  of  produced  electricity  and  heat.  Fuel  cells  behave  as  a  current 
source,  where  the  load  current  defines  the  reactant  consumption 
and,  consequently,  the  operating  conditions,  while  support  systems 
have  to  ensure  that  both  reactants  are  supplied  in  sufficient  quan¬ 
tities.  If  this  is  not  the  case,  there  is  a  high  risk  of  an  irreversible 
deterioration  [5]  of  the  cell  membrane’s  functionality,  which  will 
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eventually  lead  to  failure.  Therefore,  the  fuel-cell  control  algorithm 
has  to  adapt  to  the  varying  operating  conditions  and  at  all  times 
operate  optimally  in  terms  of  efficiency,  reliability,  and  the  health  of 
the  system,  so  as  not  to  damage  the  stack  or  shorten  its  lifetime.  An 
in-depth  discussion  on  various  approaches  and  solutions  to  stack 
control  is  given  in  [6],  while  some  other  advanced  approaches  are 
presented  in  [7,8]. 

An  important  aspect  of  introducing  fuel  cells  into  daily  use  is 
their  ability  to  respond  to  the  user’s  demands.  In  this  respect, 
fuel  cells  belong  to  the  class  of  slow-response  energy  sources 
because  of  their  complex  internal  dynamics,  related  to  mass-  and 
thermal-equilibrium  laws  inside  and  around  the  stack.  For  cer¬ 
tain  applications  this  poses  no  problems  (electromotor,  forklift, 
etc.);  however,  very  often  this  is  not  the  case.  To  increase  its  range 
of  functionality  a  fuel-cell  stack  is  usually  combined  with  a  bat¬ 
tery  (or  super-capacitor)  into  a  hybrid  power  source  [9].  Such  a 
power-generation  unit  also  consists  of  air  (oxygen)  and  hydrogen 
supply  subsystems,  temperature  and  humidity  control  subsys¬ 
tems  and  a  DC/DC  converter  [10].  During  operation  the  battery 
(or  super-capacitor)  is  used  to  sustain  the  current  peaks  that  are 
a  consequence  of  the  changes  on  the  load  side,  while  the  stack  is 
controlled  to  provide  the  right  power  at  the  right  time;  the  control 
actuator  being  a  DC/DC  converter,  which  adapts  its  output  voltage 
to  produce  the  required  currents.  A  number  of  studies  discussing 
the  control  of  such  systems  are  presented  in  [5,9-12]. 

Some  studies  follow  the  two-converter  solution,  where  they 
adapt  the  voltage  of  the  battery  and  the  stack  to  the  voltage  of 
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Nomenclature 

DC  direct  current 

FSM  finite  state  machine 

HV  heating  value 

Kp  proportional  gain  -  parameter  of  PI/PID  controller 
MOMI  magnitude  optimum  multiple  integration 
SOC  state  of  charge 

Tf  integral  time  -  parameter  of  PI/PID  controller 

Td  differential  time  -  parameter  of  PID  controller 

Greek  symbols 

X  oxygen  excess  ratio 

r]  efficiency 

cp  mass  flow 

Subscripts 

aux  auxiliary,  i.e.,  losses 

batt  battery 

comp  compressor  (blower) 

DCDC  DC/DC  converter  value 
deg  degradation 

err  error 

lim  limit/saturation  value 

meas  measured  value 

net  nett  value 

ref  reference  value 

stack  stack  value 


the  bus  [11],  while  other  approaches  [9]  (including  the  one  pre¬ 
sented  here)  assume  just  a  single  converter  between  the  stack 
and  the  other  elements.  This  partly  restricts  the  applicability  due 
to  the  somewhat  varying  output  voltage;  however,  it  has  other 
important  advantages:  simpler  control,  lower  price  and  fewer 
losses. 

In  this  work  a  control  approach  is  proposed  for  the  complete 
system.  It  is  designed  both  to  coordinate  the  operation  of  all  the 
components  and  to  optimize  the  system’s  sustainable  functioning 
in  terms  of  high  efficiency,  safety,  reliability  and  durability  [13,14]. 
The  control  is  divided  into  high-  and  low-level  control,  where  the 
high-level  algorithm  collects  the  process  data,  defines  the  current 
state  of  the  system  and  sets  references  for  the  low-level  control 
loops.  An  advantage  of  the  described  control  in  comparison  with 
others  is  its  low  complexity  and  the  fact  that  a  model  of  the  system 
is  not  required  for  tuning. 

The  paper  presents  a  simulation  system  model,  consisting  of 
a  stack,  a  DC/DC  converter,  a  battery  and  an  electric  load.  Fur¬ 
thermore,  a  control  approach  is  presented,  with  a  description  of 
the  high-  and  low-level  control  algorithms.  The  proposed  control 
strategy  is  evaluated  in  a  simulation,  where  the  results  are  also 
compared  to  another  control  strategy.  The  comparison  is  sum¬ 
marised  with  a  discussion  relating  to  the  overall  efficiency  and  the 
stack-degradation  measures. 


2.  Extended  model  of  a  FC  hybrid  power-generation  unit 

For  the  purpose  of  improving  the  operation  of  a  complete  system 
a  model  of  the  system  can  be  used.  Such  a  model  should  consist  of 
all  the  elements  in  the  energy  production-distribution  chain,  i.e., 
a  stack  model,  a  power-converter  model,  a  battery  and  a  load.  The 
scheme  of  the  complete  power  system  is  shown  in  Fig.  1. 


2.1.  Model  of  a  FC  power-generation  unit 


For  the  analysis  and  the  stack  control  design  a  model  [15]  of 
a  1.2-kW  Ballard  Nexa  power-generation  unit  [16]  was  used.  The 
Nexa  is  widely  used  in  academic  institutions  and  therefore  inter¬ 
esting  to  use  for  the  purposes  of  a  comparison.  The  model,  which 
fits  the  real  data  very  closely,  is  implemented  in  Matlab-Simulink 
and  allows  a  precise  study  of  the  system  states  and  operation. 

During  normal  operation  of  the  stack  a  (load-dependent)  suffi¬ 
cient  supply  of  reagents  (02,  FI2)  has  to  be  provided.  Flydrogen  is 
stored  in  a  high-pressure  bottle  or  metal-hydride  tanks;  its  flow 
control  is  defined  by  a  mechanical  ratio  control  that  takes  into 
account  the  required  air  flow,  and  is,  therefore,  treated  as  always 
sufficient.  The  oxygen  (from  the  air)  is  supplied  by  a  blower  that 
feeds  air  into  the  stack,  and  which  can  be  described  as  the  first  con¬ 
trol  input.  However,  for  technical  reasons,  it  has  to  be  provided  in 
excess  [5]  with  respect  to  the  chemical  equation: 

2H2  +  02  =  2H20  +  electrical  energy  +  heat  (1 ) 


Therefore,  an  oxygen-excess  ratio  X  is  defined  as: 


input  02 
reacted  02 


(2) 


The  optimal  X  value  is  dependent  on  the  stack  current  and  is  usu¬ 
ally  greater  than  2.  It  has  a  lower  critical  boundary  at  X  =  1  -  below 
this  oxygen  starvation  is  observed,  leading  to  cell  degradation  and, 
eventually,  failure.  On  the  upper  side  there  is  no  fixed  limit;  how¬ 
ever,  for  higher  values  of  X  the  compressor  requires  more  power 
than  is  gained  by  the  stack  and  so  the  output  power  decreases. 

The  air  is  also  used  to  humidify  the  membrane.  To  ensure  the 
air  brings  enough  water  it  has  to  be  either  humidified  or  taken  to  a 
temperature  where  its  relative  humidity  is  sufficient. 

The  stack  has  two  further  control  inputs.  Namely,  during  the 
operation,  besides  the  electrical  energy,  heat  is  also  produced.  The 
optimal  stack  operating  temperature  is  around  60-65  °C  [16],  and 
the  surplus  of  heat  energy  needs  to  be  taken  away  by  a  cooling  fan 
-  the  second  control  input.  Furthermore,  as  indicated  by  Eq.  (1 ),  the 
water  is  being  generated  on  the  oxygen  (cathode)  side.  This  side  is 
open,  i.e.,  there  is  a  blower  that  is  blowing  air  through  the  stack 
to  an  exhaust  where  the  water  is  released.  However,  water  is  also 
condensing  from  the  humidity  on  the  hydrogen  (anode)  side.  The 
hydrogen  side  is  a  dead  end  and  the  accumulated  water  has  to  be 
ejected  by  repeatedly  opening  the  hydrogen  purge  valve,  i.e.,  the 
third  control  input. 

The  stack  model  makes  it  possible  to  realise  the  three  described 
control  loops,  as  shown  in  Fig.  2,  i.e.,  the  water  purge  valve  control, 
the  temperature  control  and  the  oxygen  excess  ratio  control  (via  the 
compressor  voltage).  The  system  has  a  humidifier,  and  therefore  the 


Fig.  1.  Scheme  of  the  complete  power  system  -  power-generation  unit  and  load. 
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Fig.  2.  Stack  control  scheme  consisting  of  three  main  control  loops:  the  air-blower 
control  (full  line),  the  cooling/temperature  control  (dashed),  and  the  water-purge 
control  (dotted). 


humidity  can  be  assumed  to  be  at  100%.  A  comprehensive  system 
analysis  is  given  in  [17]  and  the  following  section  briefly  presents 
the  results. 

2.1.1.  Power  characteristics  of  a  fuel-cell  stack  model 

The  basic  fuel-cell  power-generation  unit  consists  of  several 
subsystems,  which  are  all,  except  for  the  stack,  direct  consumers 
of  power.  The  largest  two  power  consumers  are  the  air  blower  and 
the  cooling  fan.  Therefore,  during  a  calculation  of  the  net  output 
power  one  has  to  estimate  the  power  losses,  often  called  the  par¬ 
asitic  losses  or  the  auxiliary  power.  An  estimate  for  the  auxiliary 
power  is  proposed  in  [7],  while  in  this  work  the  following  estimate 
was  used: 

Paux  =  0.27 Ucomp  +  7.29  x  10  ^n^omp  +  1.72  x  10  4u^omp  (3) 

where  uCOmp  is  the  blower  voltage.  The  stack  unit  output  power  can 
then  be  calculated  from 


Pnet  —  Pstack  ~  Paux  (4) 

where  Pstack  is  the  stack  power,  PaUx  is  the  parasitic  loss  (auxiliary) 
power  and  Pnet  is  the  unit’s  output  power. 

To  determine  the  power  characteristics,  several  simulation 
experiments  under  various  loads  over  the  whole  operating  range 
were  made.  During  each  simulation  the  blower  voltage  was  grad¬ 
ually  decreased,  from  100%  to  30%  in  steps  of  5%.  The  results  show 
that  by  increasing  the  blower  speed  at  constant  load  (i.e.,  increas¬ 
ing  A),  the  stack  power  exponentially  (1  -k~x)  limits  its  maximum. 
However,  an  increase  of  the  blower  voltage  causes  an  increase  in 
the  parasitic  losses.  Therefore,  the  output  power  characteristics  at 
a  constant  stack  current  have  a  maximum.  The  complete  range  of 
static  characteristics  is  plotted  in  Fig.  3.  The  results  are  similar  to 
those  already  presented  in  [5]. 

2. 1.2.  Efficiency  of  a  fuel-cell  stack  model 

The  fuel-cell  stack  efficiency  is  defined  as  the  ratio  of  the  electri¬ 
cal  power  generated  by  the  stack  to  the  power  generated  by  burning 
the  same  amount  of  hydrogen  [18] 

ristacK  =  ^  Pst°V  (5) 

(Pstack.  H2hV 

where  Pstack  is  the  stack  power,  4>stackjt\2  is  the  theoretical  flow 
rate  of  hydrogen  to  the  stack,  and  HV  is  the  heating  value  of  the 
hydrogen.  In  our  case  the  hydrogen’s  higher  heating  value  (High 
HV H2  =  141 900 J  g-1 )  was  used  as  most  of  the  generated  water  is  liq¬ 
uefied.  A  small  amount  of  hydrogen  is  not  consumed  in  the  reaction 


Fig.  3.  Stack  output  power  with  respect  to  the  oxygen  excess  ratio  A,  drawn  for 
different  stack  currents.  The  points  of  maximum  output  power  are  marked  by  X. 


(due  to  the  crossover  effect  [2,3]  or  leakage  into  the  surroundings) 
and  a  modified  equation  can  be  used  in  the  efficiency  computation 

/  Pstack  $  stack  _H  2 

h  stack  ((/)stack  H2  +  (/)  stack  Joss)  ’  UV  h  stack  -f  <fistackjoss 

where  4>stackJoss  is  the  flow  of  the  actual  hydrogen  losses.  The  net 
efficiency  is  then  calculated  as: 

„  _  P net  rn\ 

n ^  ~  Supplied. H2HV  (  5 

where  Pnet  is  the  net  output  power,  and  (j) supplied  _h2  *s  suPPhed 
H2  flow  rate.  In  Fig.  4  the  calculated  characteristic  of  the  stack 
efficiency  with  respect  to  the  stack  current  is  shown,  indicating 
the  highest  efficiency  operating  point  near  Istack  =  20  A,  and  stays  at 
close  to  the  maximum  over  a  fairly  wide  range  of  current  values. 


2.2.  Simplified  power-converter  model 


There  are  two  common  types  of  DC/DC  power  converters:  the 
first  type  gives  a  lower  voltage  at  the  output  (buck  converter)  and 
the  second  gives  a  higher  output  voltage  (boost  converter).  There 
are  also  more  complex  and  expensive  converters  capable  of  both 
(buck-boost)  conversions.  In  our  case  a  buck  converter  was  used 
since  the  stack  voltage  (26-43  V)  is  always  higher  than  the  required 
battery  voltage  (24  V). 

Because  of  the  faster  dynamics,  by  several  orders,  of  the  power 
converter  compared  to  the  other  components  in  the  system,  a  sim¬ 
plified,  lumped  model  of  a  DC/DC  converter  can  be  used: 


d  = 


Uout  . 


\x  = 


Pout 


tout  Uout 
tin  Uin 


(8) 


where  d  is  the  converter  duty  cycle,  also  representing  the  out¬ 
put/input  voltage  ratio  (controlled  input),  and  (i  is  the  presumed 
constant  converter  efficiency  of  95%.  P,  U  and  I  represent  the  power, 
the  voltage  and  the  current  on  the  stack  (in)  and  the  load  (out)  side 
of  the  converter. 


Fig.  4.  Efficiency  characteristics  of  a  simulated  FC  stack  model. 
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Fig.  5.  SOC  characteristic  of  the  employed  lead-acid  battery  model  (Mat- 
lab/Simulink,  Sim  Power  Systems). 


2.3.  Battery  model 

In  our  case  a  lead-acid  battery  model  from  the  Sim  Power  Sys¬ 
tems  model  library  in  the  Matlab/Simulink  environment  was  used 
[19].  The  model  gives  the  battery  voltage  as  a  function  of  the 
input/output  currents 


Ubatt  —  E  IbattRbatt  * 


E  =  En  -K 


a-it 


+  A 


-Bit. 


t 


0 


(9) 


where  Ubatt  is  the  battery  voltage,  Ibatt  is  the  battery  current  (out), 
Rbatt  is  the  internal  resistance,  E  is  the  no-load  battery  voltage,  E0 
is  the  nominal  voltage,  Q.  is  the  battery  capacity,  I<  is  the  polariza¬ 
tion  voltage,  A  is  the  exponential  voltage  and  B  is  the  exponential 
capacity. 

The  model  also  estimates  the  battery’s  state  of  charge  (SOC), 
given  as 


SOC  =  100 


(10) 


3.  Control  design  for  the  extended  system 

Designing  a  control  scheme  for  the  entire  system  requires  a  deci¬ 
sion  on  the  hierarchy  of  the  control  objectives.  The  control  scheme 
should  comprise  three  control  loops:  (a)  charging  the  battery,  (b) 
adjusting  the  converter’s  output  voltage,  and  (c)  blowing  air  into 
the  stack.  It  has  been  shown  in  [20,21  ]  that  the  high  voltage  at  the 
stack  output,  which  translates  into  a  high  voltage  on  the  individual 
cells,  accelerates  the  growth  of  the  catalyst  particle  size,  conse¬ 
quently  reducing  its  active  surface.  At  the  same  time,  it  is  important 
to  monitor  the  lowest  cell  voltage,  which  must  not  fall  below  0.5  V 
per  cell.  In  such  a  case  it  is  possible  that  we  obtain  a  negative  voltage 
on  one  of  the  cells,  and  consequently  a  greater  voltage  on  the  neigh¬ 
bouring  cells,  which  again  leads  to  cell  deterioration.  In  addition, 
we  tend  towards  a  smooth  change  of  the  stack  load  and  operation 
in  the  range  with  the  highest  efficiency.  Regarding  the  chosen  bat¬ 
tery  type,  the  desired  control  goals  are  as  follows:  a  state  of  charge 
(SOC)  between  40%  and  80%,  charge/discharge  currents  within  the 
required  limits  and  optimization  of  the  number  of  charging  cycles. 

The  above-mentioned  requirements  lead  to  a  two-level  control 
scheme.  On  the  higher  level  there  is  an  advanced  supervisory  con¬ 
trol  algorithm,  using  the  on-line  measurements  of  both  the  battery 
SOC  and  the  actual  load  current  (J/oad),  which  calculates  a  reference 
for  the  stack  current  (I stack. ref)  with  two  goals:  maximizing  the  effi¬ 
ciency  and  ensuring  the  sustainable  operation  of  the  entire  system. 
At  the  lower  level,  two  controllers  are  used.  First,  the  DC/DC  con¬ 
verter  controller  controls  the  ratio  d,  which  affects  the  converter 
output  voltage,  uDCdc •  The  voltage  difference  ( uDCDc~Ubatt )  creates 
the  current  from  the  converter  to  the  battery,  which  further  results 
in  the  current  from  the  stack  ( htacksef  )•  The  second  lower-level  con¬ 
troller  controls  the  air  blower  so  as  to  maintain  the  appropriate 
oxygen  excess  ratio  in  the  stack  operation.  The  scheme  of  the  pro¬ 
posed  control  system  is  presented  in  Fig.  6. 

3.1.  High-level  coordination  control 

For  the  supervisory  control  algorithm  we  propose  a  control 
strategy  based  on  four  regimes  of  stack  operation 


In  our  case  the  following  nominal  values  were  used:  Q  =  10  Ah, 
E0  =  25.29  V,  Rbatt  =  0.06  Ohm,  A  =  1 .32,  and  B  =  375.  The  static  char¬ 
acteristic  of  the  battery  model  is  shown  in  Fig.  5.  In  practice,  the 
battery  capacity  is  chosen  based  on  the  stack  and  load-demand 
specifications. 


•  OFF  -  stack  shut  down. 

•  LOW  -  stack  operates  at  low  power. 

•  NORMAL  -  stack  operates  at  the  highest-efficiency  working  point. 

•  POWER  -  stack  operates  at  high  power. 


2.4.  Load  model 

The  load  may  be  generally  described  by  considering  its  purpose 
of  use.  The  latter  defines  the  energy-use  profile  and,  consequen¬ 
tially,  the  source  specification.  In  general,  we  identify  the  following 
types  of  use: 

•  Uninterruptible  power  supply  (UPS). 

o  Important  fast  start-up,  power  demand  typically  constant. 

•  Auxiliary  power  unit  (APU),  where  no  other  source  is  present. 

o  Adaptation  to  variable  load  -  variation  across  whole  operating 
range  should  not  affect  the  unit’s  life-time. 

•  Mobile  unit’s/drive  power  supply  (vehicle,  car,  forklift,  etc.). 

o  Possibility  of  a  direct  connection  of  the  stack  to  the  load  (no 
conversion  necessary). 

In  our  case  the  variable  load  was  simulated  using  a  current  gen¬ 
erator  with  instantaneous  switches,  as  the  unpredictable  on/off 
switching  of  several  power  consumers,  which  is  often  the  case  with 
APUs. 


The  switching  among  all  the  above-mentioned  regimes  is  con¬ 
ditioned  by  considering  the  instantaneous  state  of  the  battery 
(charge  =  CFIG/discharge  =  CFIG),  the  battery  SOC  and  the  load  cur¬ 
rent.  In  practice,  the  system  mostly  operates  in  the  NORMAL  regime 
at  the  optimal  working  point  (Inormai  =  20  A),  as  described  in  Section 
2.1.2.  The  LOW  mode  (J/ow  =  10  A)  is  introduced  to  prevent  switch¬ 
ing  on  and  off  too  often  in  the  case  of  low  demand.  Only  in  the  case 
of  a  high  load  current  and  a  low  SOC  value  is  the  stack  switched 
to  the  POWER  regime,  where  the  stack  current  reference  is  set 
in  accordance  with  the  actual  load  current  as  lstack.ref=m^nihoad^ 
45  A);  consequently,  with  a  lower  efficiency  and  possible  degrada¬ 
tion  effects  being  present.  In  general,  also  in  this  regime,  the  system 
tends  towards  the  lowest  possible  stack-current  values. 

The  state-transition  diagram  of  a  high-level  coordination  algo¬ 
rithm  is  depicted  as  a  finite  state  machine  (FSM)  in  Fig.  7. 
Furthermore,  considering  the  load  specifics  and  its  time  profile, 
it  is  possible  to  set  additional  time  conditions  with  respect  to  the 
switching  among  the  FSM  states,  i.e.,  the  system’s  operating  modes. 
The  actual  stack-current  reference  is  changed  following  a  ramp  pro¬ 
file.  In  this  way,  smooth  transitions  are  achieved,  ensuring  that  the 
blower  supplies  a  sufficient  amount  of  air  to  the  stack  at  all  times. 
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Fig.  6.  Two-level  control  scheme  of  the  extended  fuel-cell  system  model;  the  power  elements  have  a  gray  background;  the  full  lines  represent  current  flows;  and  the  dashed 
lines  represent  control  signals. 


3.2.  Low-level  control-DC/DC  converter  control 

The  element  that  influences  the  battery  charging  and  the  stack 
load  is  the  DC/DC  power  converter.  As  explained  in  2.2,  by  varying 
the  input/output  voltage  ratio  d,  its  output  voltage  (L^dcdc)  is  con¬ 
trolled  so  that,  consequently,  the  desired  current  Istack  is  drained 
from  the  stack.  In  our  case  the  Istack  PI  control  law  is  as  follows 

d  =  Kp(I stack  err  +  77  Jq  ^ stack  err dx)  —  -j:  f0(d  ~  d\im)dx, 

Kp  =  0.31,  Tf  =  1.148s  (11) 

where  /stacfc  err  is  the  difference  between  the  measured  value  Istack 
and  its  reference  value  ( htack-ref )•  Given  the  limited  range  of  d  val¬ 
ues  (d/im  is  represented  by  0  and  1)  the  anti-windup  protection  is 
introduced,  shown  in  the  last  part  of  the  equation.  The  controller 
was  tuned  using  the  MOMI  method  for  PI  controllers  [22]  with 
some  manual  fine-tuning.  In  this  way,  the  controller  is  tuned  with 
the  off-line-generated,  open-loop,  step-response  control  data  using 
multiple  integrations. 

3.3.  Low-level  control-air-blower  control 

A  sufficient  air  supply  is  vital  for  the  stack’s  operation.  In  order  to 
optimize  the  operation,  the  compressor  is  used  to  adjust  X  towards 
the  desired  value;  the  control  scheme  of  the  air  supply  is  shown  in 
Fig.  8. 


The  reference  X  value  of  the  air-supply  control  varies  regarding 
the  stack  current  based  on  the  maximum  power  points  shown  in 
the  characteristics  in  Fig.  3.  The  nonlinear  static  dependency  of  the 
optimum  points,  given  in  Table  1,  is  presented  in  Fig.  9. 

The  characteristic  curve  is  used  for  the  reference  X  value  calcu¬ 
lation  for  the  compressor  control  to  minimize  the  losses  and  use 
the  stack  at  the  safest  power  levels.  For  X  tracking  control  a  PID 
controller  with  the  following  control  law  was  used: 


bicomp  —  Kp(Xerr  +  Jq  Xerrdx  +  T(t  c^rr ) 

_  T’T  Jo  (“comp  —  UComp  Jim )dt, 


(12) 


where  Xerr  is  the  error  between  the  reference  and  the  estimated 
value  of  X,  and  the  controller  parameter  values  are  I<p  =  6,  Tz  =  2/3  s, 
Td  =  l/15s.  Due  to  the  limited  operating  range  of  the  compressor 
voltage  (0-100%),  the  anti-windup  protection  is  also  used  here.  The 
controller  was  tuned  using  the  MOMI  method  for  PID  controllers 
[22]  with  some  experimental  fine-tuning. 

The  nonlinearity  of  the  X  characteristic  curve  would  normally 
cause  control  problems;  therefore,  it  is  important  to  point  out  that 
the  effect  is  far  more  expressed  in  the  lower  part  of  the  operating 
range.  As  mentioned  above,  we  try  to  avoid  operating  at  low  cur¬ 
rents  as  much  as  possible  anyway,  as  it  may  lead  to  a  deterioration 
of  the  stack’s  performance  [20].  Thus,  the  controller  may  be  tuned 
to  operate  best  in  the  medium-high  operating  range. 


Fig.  7.  Regimes  and  conditions  of  the  proposed  (Finite  State  Machine)  high-level  controller. 
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Fig.  8.  Air-supply  control  scheme. 


Table  1 

Static  dependency  of  optimum  k  values  according  to  the  stack  current  ( IstaCk  )• 


htack  (^) 

5 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

k 

3.522 

2.662 

2.375 

2.230 

2.143 

2.084 

2.042 

2.010 

1.985 

1.965 

1.947 

Fig.  9.  Optimum  k  values  for  the  entire  stack  operating  range. 


Fig.  10  shows  the  time  plot  for  the  stack  signals  for  the  vary¬ 
ing  load  profile  in  the  upper  diagram,  for  a  safe  value  Xref=3  (gray 
line)  and  for  a  power  optimized  Xref=Xopt  (black  line),  where  a 
linear  interpolation  between  the  Xopt  points  from  Table  1  is  used. 
Using  Xopt  values  in  a  control  algorithm,  a  slightly  higher  net  out¬ 
put  power  in  the  second  diagram  and  a  higher  net  efficiency  in 
the  third  diagram  can  be  observed,  especially  in  the  higher-power 
operating  range.  However,  in  both  cases  the  key  control  problem 
is  providing  a  sufficient  quantity  of  air  at  the  transients,  where 
the  consumption/load  is  suddenly  increased,  and  consequently  X 
may  temporarily  fall  below  the  critical  value  of  1,  as  happens  at 
1500s. 


time  [s] 


Fig.  10.  Stack  control  with  a  variable  load  profile  for  reference  k  =  3  (gray)  and  k  =  kopt  (black).  Diagrams  top-down:  (1 )  stack  current,  (2)  net  output  power,  (3)  net  efficiency 
and,  (4)  oxygen  excess  ratio  k. 
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Table  2 

Operating  conditions’  impact  on  stack  degradation  for  a  1.2-kW  Nexa  power  unit  produced  by  Ballard. 
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Operation 

Conditions 

Degradation  Rate 

Period 

Losses 

Steady-state 

Full  power  (46  Ane(:) 

-25  mWh-1 

1500h 

-37  W 

-0.54  mV  h"1 

-0.8  V 

Part  load  (<35  Anet)a 

-OmWh-1 

1500h 

OW 

-OmVh-1 

0.0  V 

Dynamic  loads 

Idle  to  full  power 

-OmWh-1 

1500h 

OW 

-OmVh-1 

0.0  V 

On/off 

o 

i 

CTl 

> 

-52  mW  cycle-1 

500  cycles 

-26  W 

-1.10  mV  cycle  h  1 

-0.56V 

a  We  considered  that  the  losses  increase  linearly  from  0  to  25  mWh-1  (i.e.,  0.54mVh_1)  as  the  current  increases  from  35  A  to  a  maximum  of  46  A. 


4.  Criteria  for  evaluation  of  the  proposed  control 

A  useful  comparison  of  energy  sources  can  only  be  based  on 
measurable  quantities.  Below  we  present  the  two  criteria  that  will 
help  us  in  assessing  the  overall  quality  of  the  control.  These  criteria 
are  the  stack  efficiency  and  stack  life-time  degradation  impact. 


4.1.  Stack  efficiency 


One  of  the  key  parameters  for  the  stack-operation  evaluation  is 
efficiency.  In  this  area  the  fuel  cells,  as  an  energy  source,  compete 
with  others.  To  evaluate  the  efficiency,  we  used  a  simple  measure 
-  the  integral  of  the  ratio  of  the  generated  output  electrical  power 
(less  the  auxiliary  power)  and  the  power  of  burning  the  hydrogen 
fed  to  the  stack  (see  Eq.  (7)): 


Coverall 

Jo 


Tlnetdt 


(13) 


4.2.  Stack  life-time  degradation 

To  assess  the  influence  of  the  operating  conditions  on  the  deteri¬ 
oration  of  the  stack’s  characteristics,  we  used  the  information  from 
documentation  [16]  that  the  manufacturer  provides  for  the  1 .2-kW 
Nexa  power  unit.  They  have  tested  specific  conditions  for  1500h 
and  500  off-to-full  power  cycles,  for  which  the  data  are  shown  in 
Table  2. 


Table  3 

Degradation  influence  of  the  oxygen-starvation  effect. 


Operation 

Conditions 

Degradation  rate 

Steady-state 

A,  <  1.0 

A  =  1.0a 

-OmWs-1 

A  =  0.0 

-lOrnWs-1 

a  We  considered  that  the  losses  increase  linearly  from  0  to  lOrnWs-1  as  the  A 
decreases  from  1  towards  0. 


As  an  additional  degradation  measure,  the  degradation  rate  for 
stack  operation  below  the  oxygen  excess  ratio  X  =  1  was  introduced, 
given  in  Table  3.  Cell  starvation  is  a  very  rigorous  degradation  phe¬ 
nomenon  and  therefore  it  is  heavily  penalized.  This  data  is  not 
given  in  the  Nexa  documentation  as  the  embedded  controller  and 
DC/DC  converter  provided  in  the  package  work  to  avoid  such  con¬ 
ditions.  The  influence  escalates  linearly  with  the  distance  from 
value  1. 

As  an  indicator  in  the  control  strategy  comparison  we  used 
the  cumulative  sum  of  the  performance  losses  due  to  operation 
in  degenerating  conditions  over  the  simulation  time. 


Pdeg  ~ 


Ustack  ~  35)  25  x  10~3 
46  -  35  3600 


+  max  (  fj — 10x10 


'-3,o) 


dr  +  n  x  52  x  10 


-3 


(14) 


where  Pdeg  is  the  cumulative  power  fall-off  as  a  consequence  of  the 
degradation  and  n  is  the  number  of  on/off  cycles. 


Fig.  11.  Simulation  scheme  of  the  extended  system  model  and  control  in  the  Matlab/Simulink  environment. 
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Fig.  12.  Extended  system-control  time  plot;  a  constant  DC/DC  converter  transmission  ratio  dDcDc  =  0.75  is  used.  Signals  top-down:  (1)  load  current  (dashed  black),  stack 
current  (black),  DC/DC  out  current  (gray);  (2)  stack  voltage  (black)  and  DC/DC  out  voltage;  (3)  DC/DC  transmission  ratio  d  (gray)  and  battery  SOC  (black);  (4)  efficiency  (black) 
and  degradation  (gray);  (5)  reference  (black)  and  actual  (gray)  oxygen  excess  ratio  values. 


5.  Control  strategy  analysis  under  a  simulated  load  profile 

The  proposed  approach  of  the  extended  FC  system  control  was 
compared  to  another  approach  in  a  simulation  where  a  constant 
value  of  the  DC/DC  power  converter  voltage  ratio  d  is  assumed.  Its 
value  is  defined  by  the  average  ratio  of  the  stack-to-battery  voltage. 

The  simulation  scheme  of  the  extended  model  (presented  in 
Section  3)  and  the  proposed  control  system  using  Matlab/Simulink 
is  shown  in  Fig.  11.  In  the  simulation  a  smaller-capacity  battery 
(Q  =  1 0  A  h)  was  used,  so  that  in  a  reasonable  time  some  interesting 
responses  and  transients  could  be  observed. 

In  Figs.  12  and  13  the  simulated  essential  signals  of  the  evalu¬ 
ated  control  strategies  are  presented.  The  experimental  load  profile 
is  plotted  with  a  dashed  line  in  the  upper  diagram.  In  Fig.  12  a  sim¬ 
ulation  using  a  constant  DC/DC  conversion  ratio  (d  =  0.75)  is  used. 
It  can  be  observed  that  there  is  no  direct  control  over  the  bat¬ 


tery  SOC  (Diagram  3),  and  its  value  falls  down  to  a  critical  level 
(12.6%)  during  a  very  high  load;  similarly,  the  overcharging  prob¬ 
lem  arises  with  small  loads.  In  addition  the  load  variations  (steps, 
transients)  transform  directly  to  stack-current  spikes,  resulting  in 
A  (Diagram  5)  falling  below  the  value  of  1  on  certain  occasions, 
consequently  leading  to  stack  degradation.  On  the  other  hand, 
the  operating  point  is  kept  close  to  the  optimal  efficiency  at  all 
times. 

In  Fig.  13  the  simulation  using  a  control  approach  that  is 
described  in  this  work  is  presented.  One  can  see  that  by  using  the 
DC/DC  converter  and  battery  it  is  possible  to  compensate  for  the 
current  peaks  and  spikes  caused  by  the  load  changes  and  achieve 
smoother,  ramp-shaped  transitions  of  the  stack  current.  The  stack  is 
controlled  to  operate  in  optimal  conditions  whenever  possible,  only 
increasing  the  stack  power  in  the  case  of  large  loads  and  decreasing 
it  when  the  battery  is  highly  charged  and  the  load  current  is  small. 
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Fig.  13.  Extended  system-control  time  plot  of  the  strategy  that  is  presented  in  this  work.  Signals  in  the  top  Diagram  (1 )  load  current  (dashed  black),  stack-current  reference 
(dashed  gray)  and  actual  (black)  value,  DC/DC  out  current  (gray).  Signals  in  other  diagrams  as  in  Fig.  12. 


The  reader  should  note  that  the  SOC  values  after  the  experiment 
are  not  equal  for  the  presented  control  approaches.  This  is  because 
the  load  profile  was  considered  as  an  excitation  signal  only  and  the 
study  focuses  on  the  stack-operation  parameters  and  the  quality  of 
the  energy  provided.  The  criteria  values  for  the  control  strategies 
are  given  in  Table  4. 

6.  Discussion 

From  these  simulation  experiments  some  interesting  conclu¬ 
sions  can  be  drawn.  In  terms  of  efficiency  only,  there  is  an  advantage 
with  the  simpler  approach  using  a  constant  d  ratio  over  the  con¬ 
trol  approach  presented  in  this  work.  Using  a  constant  d,  a  subtle 
closed-loop  control  function  over  the  battery  SOC  is  embedded.  On 
average  it  causes  the  stack  current  to  tend  towards  the  appropriate 
constant  value  (see  Fig.  12),  provided  d  is  set  correctly,  and  the  load 
profile  has  a  constant  long-term  average.  In  the  presented  case  the 


choice  of  d  =  0.75  resulted  in  a  lower  final  SOC  value  and  a  lower 
average  output/load  voltage  compared  to  the  supervisory  control; 
by  increasing  the  d  ratio  the  efficiency  would  decrease.  However, 
using  constant  d,  the  charging  depends  mainly  on  the  instantaneous 
load  and  the  SOC,  so  the  desired  charge/recharge  cycles  and  SOC 
limits  cannot  be  ensured. 

This  shows  that  by  designing  a  sensible  control  function  a  bal¬ 
ance  between  both  approaches  can  be  achieved,  which  would  allow 
a  further  optimized  operation  to  be  achieved  without  sacrificing 
the  safety  of  the  components  during  operation.  Another  point  to 
note  is  that  including  knowledge  of  the  future  load  profile  would 


Table  4 

Comparison  of  indicators  for  the  tested  control  strategies. 


d  =  0.75 

FSM  supervision  control 

Net  stack  efficiency 

0.4099 

0.3804 

Stack  degradation  (mW) 

31.1286 

7.0773 

9428 


B.  Pregelj  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  9419-9428 


make  it  possible  to  plan  the  operation  closer  to  optimal  conditions 
and  bring  additional  gains  in  terms  of  efficiency  maximisation  and 
component-degradation  reduction. 

7.  Conclusion 

We  present  a  model  of  a  fuel-cell-based  power-generation 
unit,  consisting  of  the  stack,  a  battery  and  a  power  converter.  An 
advanced,  two-level  control  approach  is  presented,  consisting  of  a 
finite  state  machine  for  the  high-level  control  and  PID  controllers 
for  the  low-level  loops,  which  provides  a  good  balance  between 
optimal  and  reliable  control,  as  well  as  a  low  computational 
complexity.  This  approach  to  the  design  of  FC  system  control 
is  compared  to  a  simpler  strategy  and  evaluated  with  efficiency 
and  degradation  criteria.  The  proposed  control  achieves  a  slightly 
smaller  efficiency,  but  maintains  the  stack  and  the  battery  within 
the  safe  operating  range.  In  future  investigations,  a  more  advanced 
control  algorithm  with  a  relatively  simple  model  of  the  complete 
system  will  be  developed  and  tested  on  a  real  system. 
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